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ABSTRACT 



We present new multi-epoch near-infrared (NIR) and optical high-angular 
images of the V 773 Tau pre-main sequence triple system, a weak-line T Tauri 
(WTTS) system in which the presence of an evolved, "fossil" protoplanetary 
disk has been inferred on the basis of a significant infrared (IR) excess. Our 
images reveal a fourth object bound to the system, V 773 Tau D. While it is 



much fainter than all other components at 2 /im, it is the brightest source in the 
system at 4.7/xm. We also present medium-resolution K band adaptive optics 
spectroscopy of this object, which is featureless with the exception of a weak Br7 
O ' emission line. Based on this spectrum and on the spectral energy distribution 

(SED) of the system, we show that V 773 Tau D is another member of the small 
class of "infrared companions" (IRCs) to T Tauri stars (TTS). It is the least 
luminous, and probably the least massive, component of the system, as opposed 
^ ■ to most other IRCs, which suggests that numerous low-luminosity IRCs such as 

V 773 Tau D may still remain to be discovered. Furthermore, it is the source of 
the strong IR excess in the system. We therefore reject the interpretation of this 
excess as the signature of a fossil (or "passive") disk and further suggest that 
these systems may be much less frequent than previously thought. 

We further show that V 773 Tau C is a variable classical TTS (CTTS) and that 
its motion provides a well constrained orbital model. We show that V 773 Tau D 
can be dynamically stable within this quadruple system if its orbit is highly 
inclined. Finally, V 773 Tau is the first multiple system to display such a variety 
of evolutionary states (WTTS, CTTS, IRC), which may be the consequence of 
the strong star-star interactions in this compact quadruple system. 



Subject headings: stars: individual (V 773 Tau) — binaries: close — stars: pre- 
main sequence — planetary systems: protoplanetary disks 
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1. Introduction 

V 773 Tau (HD 283447, HBC 367) is a remarkable ~ 5 Myr-old, ~1.5M (White & Ghez 
2001) pre-main sequence object located 148 pc away (Lestrade et al. 1999) in the Taurus 
star-forming region. Its SED displays two broad maxima at ~ 1 jim. and ~ 20 /im (cf, e.g., 
data compilation by Feigelson et al. 1994). The first one corresponds roughly to what is 
expected for a stellar photosphere, while the longer wavelength component is reminiscent of 
the spectral energy distribution (SED) of accreting pre-main sequence stars, i.e. classical 
TTauri stars (CTTS). This infrared (IR) excess in CTTS is generally attributed to a dusty 
circumstellar disk heated by reprocessing of starlight and by energy released through the 
accretion process. While the IR excess in V 773 Tau suggests a disk model for this system, 
the weak hydrogen emission lines in its spectrum (e.g., Cabrit et al. 1990) implies that no 
accretion is ongoing on the stars, consistent with a weak-line TTauri star (WTTS), although 
the presence of forbidden emission lines in its spectrum imply mass loss. This has been 
considered as substantial evidence for the presence of a "passive disk" around the stars, i.e., 
a non- accreting dusty structure that only reprocesses the stellar light. Furthermore, Phillips 
et al. (1996) first pointed out that, because of the dip around 2 /im in the SED of V 773 Tau, 
this object presents the characteristics of a "fossil disk", as defined by Strom (1995), in 
which the inner parts of the disk have been cleared out of small dust grains possibly through 
the formation of planetesimals. In other words, this system has been proposed to be one of 
the few "transition objects" between CTTS and WTTS, observed during the phase of disk 
dissipation. This duality WTTS/IR excess is present in only a few objects (e.g., Skrutskie 
et al. 1990; Stassun et al. 2001), which suggests that this is a short-lived phase. 

While the analysis of the SED and the emission lines have been based on low angular 
resolution data sets, high angular and spectral resolution studies have resolved this object 
into a triple system. First, a tight (~ O'.'l 15 AU) stellar companion was discovered through 
speckle imaging (Ghez, Neugebauer & Matthews 1993; Leinert et al. 1993). The primary 
component of this pair later appeared to be itself double, a 51 day-period (~ 0.3 AU semi- 
major axis; Welty 1995) double-lined spectroscopic binary (SB). All three components are 
believed to be WTTS (White & Ghez 2001). Interactions between the various components 
are likely to be important in this tight multiple star system. Indeed, the variable non- 
thermal radio emission observed in this system (O'Neal et al. 1990; Dutrey et al. 1996) 
appears to be phase- locked with the SB (Massi, Menten & Neidhofer 2002). This emission 
therefore appears to be chromospheric with enhanced activity occuring when the two stellar 
chromospheres partially overlap at periastron. The properties of the putative fossil disk 
could also be affected by the multiplicity of this object. A fit of a disk emission model to 
the IR excess places the inner radius of the optically thick disk at about ~ 0.2 AU. This is 
well inside the orbit of the SB, raising doubts about the disk model. 
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In this paper, we present new near-infrared (NIR) and optical high-angular resolution 
imaging of V773Tau, which unambiguously reveal a fourth component. The observed 1.6- 
4.7/xm colors and the broadband 2 fim spectrum of this object show that it is an infrared 
comapnion (IRC) and that it is in fact the source of the IR excess in this system. Throughout 
this paper, we will refer to the SB as V773TauAB, the long-known speckle companion as 
V773TauC and the new companion identified in this study as V773TauD. We present 
the various datasets used in this paper in §2. The orbital motion of V773TauC and the 
photometric properties of all components of the system are analyzed in §3. In §4, we 
determine the status of each component and study the importance of star-star interactions 
in the system. Finally, we summarize our results in §5. 



2. Observations and data analysis 

2.1. Speckle imaging 

V773Tau has been regularly targeted by speckle interferometry surveys because, as 
a ~ O'.'l binary, it is an excellent candidate for studies of orbital motion. A number of 
measurements (see Table 1) have already been reported (Leinert et al. 1993; Ghez et al. 
1993, 1995, 1997; Woitas, K6hler & Leinert 2001, Tamazian et al. 2002). 

We present here several new NIR speckle observations of the system. These data were 
obtained at 2.2 ^m (K band) on 1996 Oct 24, 1997 Nov 19, 1998 Nov 04, 1999 Nov 20 
and 2000 Dec 08 at the 200-inch telescope at Palomar Observatory and on 2001 Dec 08 at 
the 10 m Keck I telescope. We used a 64x64 subsection of a 256x256 SRBC InSb aray at 
Palomar and the facility instrument Near-IR Camera (Matthews & Soifer 1994; Matthews 
et al. 1996) at Keck. These instruments offer pixel scales of 0'.'0337±0'. / 0005/pixel and 
0'. / 0203±0'/0003/pixel respectively, adequately sampling the diffraction limit in the K band 
(~ 0'.'09 and ~ 0'.'05), and have known orientations to an accuracy of ±1°. Stacks of several 
hundred short integration (t ~ 0.1s) exposures were obtained, immediately followed by a 
similar stack of exposures on a calibration point source. Standard speckle data reduction 
routines were applied to the data; we refer the reader to Ghez et al. (1993) and Patience et 
al. (1998) for more details. Examples of power spectra at various epochs are presented in 
Figure 1. 

In the latter four epochs, the calibrated power spectra of the system reveal two distinct 
sets of fringes, corresponding to the V773TauAB-C and V773TauABC-D pairs 1 . In those 



1 The second set is in fact the sum of the V773TauAB-D and V773TauC-D almost parallel sets of 
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cases, we first fit the better-defined V773TauAB-C pair, then divide the observed power 
spectrum of the system by a model of this binary and finally fit a second binary system. 
The information regarding the V 773 Tau AB-D pair is then easily recovered via the known 
separation and flux ratio of V773TauAB-C. A regular binary system model is fit to the 
other datasets, which do not show sufficient evidence for a third resolved component 2 . 

We also present here the results of a visible wavelength speckle interferometry exper- 
iment performed on 1991 Sep 25 to 29 at the 200-inch telescope at Palomar Observatory 
with the same setup as that described by Gorham et al. (1992). The filters used, Fq 2 o, -^656; 
Fgoo and Fsso, are all 100 A wide and are centered on 6200, 6560, 8000 and 8500 A respec- 
tively. Data reduction and analysis routines are also similar to that described by Gorham 
et al. (1992), in which autocorrelation functions (ACFs) are generated. While ACFs are 
related to the power spectrum analyzed for the IR data, they are more robustly generated 
in photon-counting experiments. Only two stars are identified in our data, which are readily 
associated with V773TauAB and C. 



2.2. Adaptive optics imaging 

2.2.1. H & K band imaging 

In addition to this wealth of speckle measurements, we also obtained H band 3 and K 
band adaptive optics (AO) images of the V 773 Tau system with the AO system of the 10 m 
Keck II telescope (Wizinowich et al. 2000) and the facility NIR spectrometer NIRSPEC 
(McLean et al. 2000). The data were obtained on 2000 Nov 20 with SCAM, the slit- 
viewing camera inside NIRSPEC. With each filter, we obtained 4 sets of 100 individual 0.1 s 
exposures, with the object centered on a different part of the detector in each set. The 



fringes, which we treat as a single fringe pattern corresponding to V 773 Tau D and the photocenter of 
V 773 Tau ABC. 

2 In the 1997 dataset, the V773TauAB-C is not clearly resolved although the amplitude of the power 
spectrum decreases towards the outer edge of the power spectrum along PA^ 40° (see Figure 1). Since we 
do not detect the first minimum of the visibility curve along this direction, we can only set an upper limit on 
the binary separation (< X/2D « 45 mas in the K band with the Palomar 200- inch telescope). In subsequent 
epochs, this first minimum is reached and the binary parameters can be confidently derived. 

3 Thc filter we used, "NIRSPEC-3" , is a customized filter specific to the instrument. Its central wavelength 
(A c = 1.625 /im) and bandwith (AA = 0.36 /im) are close enough to the standard broadband H filter 
(A c = 1.630 /im, AA = 0.31 /im) that we do not expect magnitude differences larger than a few hundredths. 
Throughout our paper, we therefore refer to this filter as the H filter. 
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individual exposure time was short enough to avoid saturating the stellar cores. V 773 Tau 
is bright enough in the visible (V ~ 10.7, R ~ 9.8) and is unresolved in the wavefront sensor 
camera so it could be used as a natural guide star for the AO system. We achieved Strehl 
ratios of 12% and 15% respectively at H and K, which is enough to provide diffraction- 
limited cores. The resolution of the images, measured as the FWHM of a Gaussian fit to the 
radial profile of the brightest component, V773TauAB, are about 40 and 50 mas at H and 
K, respectively. The final images are presented in Figure 2 (upper row). 

Data reduction included sky subtraction, flat-fielding, bad pixel removal and shift-and- 
adding the individual images and was performed in IRAF 4 . Astrometry and relative pho- 
tometry are performed through two methods. First, we estimate the centroid locations for 
each component in the final image. To estimate uncertainties, we perform the same analysis 
on four subsets of the data and calculate their standard deviation. Relative photometry is 
obtained by summing fluxes over 2-pixel radius apertures around each component; uncertain- 
ties are similarly obtained through subsets of the data. To avoid potential star-to-star flux 
contamination, we also fit point spread functions (PSFs) to the multiple system. Because 
of the limited field of view of SCAM, no other star is available to be used as a contempora- 
neous PSF. We thus use both a theoretical image which is basically the Fourier transform 
of the entrance pupil of the telescope and an empirical PSF with a somewhat better image 
quality, which we obtained on Nov 19 by imaging TTauN (see Duchene, Ghez & McCabe 
2002 for details). All three techniques agree within 3a of each other and all trends are 
naturally explained by the different imperfections in the PSFs. The final astrometric and 
photometric results presented in Table 1 are the averages of these three techniques (one direct 
measurement and two fits). 



2.2.2. L' & M' band imaging 

We also imaged vt at 3.8 fim and 4.7/im (L' and M' broadband filters, which are 0.7/xm 
and 0.24 /xm wide respectively) with the new facility instrument NIRC2 (Matthews et al., 
in prep.) behind the AO system on the Keck II telescope. The data were obtained on 
2002 December 13 under excellent seeing conditions; V 773 Tau itself was used as a natural 
AO guide star. Strehl ratios on order of 60-75 % were measured at L' on similarly bright 
single stars throughout the night. A circular cold mask equivalent to a circular 9 m-diameter 
entrance pupil was used to minimize the background emission, a critical factor in the thermal 



4 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under contract to the National Science Foundation. 
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IR range. The resulting diffraction-limited images have Gaussian FWHM of 90 and 110 mas 
at V and M', respectively. 

Several hundred short exposures (~ 0.05 s each) were obtained at four different locations 
on the detector for a total integration time of about lmin with each filter. The datasets 
were subtracted from each other to remove the thermal background, flat-fielded, corrected 
for bad pixels and shift-and-added to produce the final images presented in Figure 2 (lower 
row). Astrometric and photometric parameters of the multiple systems were obtained using 
the same methods as described above for the H and K images. The PSF fitting process is 
performed using HD 1160 (a single star) at V and T Tau N at V and M' . Both were observed 
on the same night with the same setup. All methods yield compatible estimates within ~ 2 a 
for both the location of the centroids and the flux ratios. 



2.3. Adaptive optics spectroscopy 

On 2000 Nov 20, we also obtained long-slit K band spectroscopy of two components of 
the system, by aligning the 07036-wide slit of NIRSPEC along the V 773 Tau AB-V 773 Tau D 
position angle. Four spectra were obtained at different locations behind the slit, with a total 
integration time of 10 minutes. The instrumental setup is identical to that used in our 
previous study of the tight binary TTauS reported in Duchene et al. (2002), which also 
provides more details on the data reduction process. The spectra cover the 2.0-2.4 //m 
range, with a 2-pixel resolution of R ~ 2500. Immediately after V 773 Tau, we observed the 
nearby stars HD 27311 (AO) and HD 27560 (GO) to correct for telluric absorption features. 
The spectra were first sky-subtracted, corrected for distortion effects, flat-fielded and bad 
pixel-corrected. 

To extract the individual spectra, we performed a fit to the binary profile along the 
spatial axis with a semi-empirical PSF. The latter is obtained by unfolding the Northern half 
of the profile of V 773 Tau AB, i.e., assuming a symmetric profile. The same profile was fitted 
to both components simultaneously. We estimate that the contamination of the spectrum 
of V 773 Tau D by that of V 773 Tau AB is less than 5 % of the former at any wavelength. 
The continuum around Br7 in the spectra of V 773 Tau was not perfectly corrected with 
HD 27560, either because its spectrum is not close enough to the solar spectrum or because 
the telluric features changed between our observations. We thus divided both spectra by a 
third-order polynomial function fit to a 40 nm-wide window centered on Br7 and containing 
no other photospheric feature. The resulting spectrum of V 773 Tau AB in this window 
(except for about ten pixels around the hydrogen line itself) was used to correct the spectrum 
of V 773 Tau D of any remaining telluric feature. This procedure does not modify the shape 
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or strength of the Br7 line in either spectrum but it improves significantly the continuum 
around that line. Finally, all spectra were divided by a fifth-order polynomial function fit to 
their respective continua. The output spectra are presented in Figure 3. The signal-to-noise 
ratio is estimated to be on the order of 100 per pixel in the final spectra. 

3. Results 

3.1. A new stable quadruple system in Taurus 

As can be seen in Figure 2, the system is clearly resolved in all our AO images into three 
components at a spatial resolution of O'.'05-O'.'l. The brightest component at 1.6-3.8 /xm also 
dominates the flux in the visible and is associated with V773TauAB, the SB discovered by 
Welty (1995). The latter cannot be resolved in our data; even at apoastron its projected sep- 
aration on the sky is 0'.'0029. The three visual components are also simultaneously detected 
in the latest epochs of speckle imaging, from 1998 onward (see Figure 1). We attempted to 
re-analyze earlier datasets, but the fainter component did not show up because of the lower 
quality of earlier IR detectors, variability of the source and/or a less favorable configuration 
of the system. Our multiwavelength dataset unambiguously reveals that V773TauAB has 
two distinct apparent companions. 

By combining our new astrometric results with previous measurements, we can study 
the relative motion of these stars. We choose to use the rest frame of V773TauAB. We 
exclude the 1991 measurement of Leinert et al. (1993) and the 1993 J band measurement of 
Woitas et al. (2001) from our analysis. The former results are obtained from deprojecting 
the visibilities measured along two orthogonal axes while the latter are obtained with an 
instrument setup that does not provide Nyquist sampling of the PSF thus creating aliasing. 
Furthermore, to minimize the number of instrumental set-ups used in this analysis without 
reducing the time coverage of the orbit, we do not include the measurements presented in 
Tamazian et al. (2002). The locations of the companions as a function of time are illustrated 
in Figure 4. The motion of both visual companions to V773TauAB are almost linear over 
the truncated 4 year time baseline 1998-2002, with velocities of 18.0±0.4mas/yr at position 
angle 125?1±2?1 and 13.7±1.1 mas/yr along position angle 120°±8° for V773TauC and 
V 773 Tau D respectively. While the proper motion of the system is of similar amplitude, its 
position angle is different (~ 152°, Frink et al. 1997), and it is unlikely that V 773 Tau D 
is a foreground/background object unrelated to the other components of V 773 Tau. This 
is further reinforced by the very small separation of the V 773 Tau AB-D pair. The surface 
density of K < 9.0 objects (since K AB ~ 7.0, Ghez, White & Simon 1997, we estimate that 
K D ~ 9.0) around V 773 Tau in the 2 Micron All-Sky Survey is about 2.5 x 10 -6 per square 
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arcsecond. The probability of finding such a foreground/background object within 0'/25 of 
V 773 Tau AB is only « 5 x f CT 7 . We therefore conclude that V 773 Tau D really is physically 
associated with V 773 Tau, which is therefore a quadruple system. 

With measurements regularly sampling the last twelve years, the orbit of V 773 Tau C 
is well defined. We therefore fit a Keplerian orbit to our astrometric dataset to constrain 
the orbital parameters, in particular the total system mass. We neglect here the influence 
of V773 TauD on the V773 Tau AB-C system. We assumed that the orbit is closed (e < 1), 
fixed the distance to D — 148 pc, and explored the 7-dimension parameter space (P, a, e, 
i, u, Q, To) using a x 2 minimization routine. We used 10000 different starting points that 
were randomly selected within the allowed range for each orbital parameter to make sure 
that we converge on the absolute minimum of x 2 ■ The resulting best fit orbit (x 2 = 1.94) is 
shown in Figure 4 and the corresponding orbital parameters are presented in Table 2. The 
estimated inclination is very similar to the orbit of the SB (Welty 1995), suggesting a coplanar 
triple system. The derived periastron and apoastron distances are 14 ± 1 AU and 26 ± 4 AU 
respectively. Until relative radial velocities within the system are measured, an ambiguity 
on the inclination cannot be solved and, independently of this, the total system mass is only 
determined to within a distance- related scaling factor: Mabc — (3.7±0.7) x (D/1A8 pc) 3 M & . 
Therefore, the uncertainty on the distance to the system (Lestrade et al. 1999) introduces an 
additional 10 % uncertainty on the system mass. Surprizingly, the uncertainties we estimate 
on the orbital parameters are not smaller than those derived by Tamazian et al. (2002) in 
spite of our extended coverage of the orbit. Some parameters differ by up to 7-9 a (P or b, 
a) between the two orbital solutions, suggesting that Tamazian et al. underestimated their 
uncertainties. 

The subsystem V 773 Tau ABC is a hierarchical triple system with a ratio of orbital 
periods of Pab-c/Pa-b ^ 300; it is therefore expected to be dynamically stable over long 
timescales. On the other hand, the projected separation of V 773 TauD is only slightly 
larger than the apoastron distance of V 773 Tau C, which raises the question of the stability 
of its orbit. Since V773TauAB has such a compact orbit, we can consider it a single point 
mass. We apply the following criterion, derived by Eggleton & Kiseleva (1995) for a non- 
coplanar triple system to V773 Tau AB-C-D: for a given combination of the mass ratios in 
a hierarchical system, there is a critical value of the ratio of the periastron distance of the 
outer pair to the apoastron separation of the inner binary beyond which the system can 
be considered stable for at least several hundred times the outer period. If this criterion 
is not fulfilled, there are times at which the three star-to-star distances in the system are 
commensurable, which is a highly unstable configuration. Applying this criterion with q in = 
M A b/Mc « 3.7 and q out = M A bc/M d > 3.1 (see §4.1), we find that the orbit of V773TauD 
is stable if its periastron distance is larger than about 115 AU. Since its projected separation 
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is about 30 AU, this suggests that its orbit is highly inclined to the line of sight (i > 75°, 
similar to the ~ 67° inclination of both V 773 Tau AB and V 773 Tau C). If V 773 Tau D were 
on a circular orbit with such a radius, its predicted orbital velocity would be ~ 6km/s. 
This is very similar to its observed projected velocity with respect to the center of mass 
of V773TauAB-C, Vb = 6.9 ± 0.8km/s as measured from a linear fit to the motion of 

V 773 Tau D over the 1998-2002 period. Our astrometric dataset is therefore consistent with 

V 773 Tau D being on a stable, inclined orbit around the inner triple system. 

3.2. SED and variability of the components of V 773 Tau 

In an effort to understand the nature of V 773 Tau D, we first construct the SED of the 
system from 4000 A up to 60 pm. Unresolved photometry for the quadruple system is taken 
from White & Ghez (2001) in the visible, Rydgren & Vrba (1981) in the 1-4 /xm range and 
Prusti et al. (1992) for IRAS measurements 5 . Flux ratios within the system are from White 
& Ghez (2001) in the visible and from this study (our AO measurements) from 1.6 //m to 
3.8 jum. No absolute photometric measurement of V 773 Tau is available at 4.7/xm, so we 
did not use the results of our M' image in this analysis. Finally, we used the R band flux 
ratio of the SB from Welty (1995) to separate the contributions of V 773 Tau A and B. To 
correct for the foreground interstellar extinction, all measurements were dereddened with 
Ay = 1.39 mag (Welty 1995, White & Ghez 2001) using a standard extinction law (Rieke 
& Lebofsky 1985). We assumed that none of the components of V 773 Tau was affected by 
circumstellar extinction. The resulting SED is presented in Figure 5. 

We fit photospheric SEDs (as tabulated by Kenyon & Hartmann 1995) to the two 
components of the SB using the VRI photometry since this wavelength range is the least 
likely to be affected by accretion- induced excesses that can be present in TTS. Initially, the 
spectral types are considered to be those derived in the literature, K2 and K5 (Welty 1995), 
but a marginally better fit is obtained with spectral types K2 and K7 for V 773 Tau A and 
B respectively; the associated stellar luminosities are 2.2 and 1.4 L , with an uncertainty of 
about 0.1 Lq. White & Ghez (2001) obtained a similar value for the luminosity of V 773 Tau A 
although they estimated it based on the / band photometry only. We do not find evidence 
for any excess above photospheric levels up to 3.8 /xm for this pair. The unresolved spectrum 
of this pair displays many photospheric features, including a weak Br7 absorption line (see 



5 Thc IRAS observations are blended with the CTTS FM Tau, which is located 40" away. Prusti et al. 
(1992) were able to disentangle the contributions of both sources at 12 and 25 fim. At 60 fim, they assumed 
the same flux ratio as observed at 25 /xm. This extraction process results in large error bars (15 % and 20 % 
at 25 and 60 /im respectively) . 
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Figure3). From a comparison of this spectrum with published spectra libraries (Kleinmann 
& Hall 1986; Wallace & Hinkle 1997), we estimate that the strength of the various features 
is compatible with a spectral type in the range K0-K5IV, consistent with the spectral types 
adopted here. 

As already pointed out by White & Ghez (2001), the SED of V773TauC can be fit 
with a pure photosphere from U to K. However, a ~ 3<r excess is then present at 3.8 /zm, 
which is further reinforced by the fact that V 773 Tau C becomes brighter than V 773 Tau AB 
at M'. The nature of V 773 Tau C is revisited in §4.1, where we argue that it is in fact a 
CTTS. For now, we only fit a photosphere to the visible measurements, and find a best fit 
for a spectral type MO. 5 with a 1 subclass uncertainty and a 1.5 ± 0.2 L Q stellar luminosity, 
also in agreement with those obtained by White & Ghez (2001). This results in a ~ 0.2 mag 
excesses at both H and K, a ~ 0.7 mag excess at V and a K — V R3 0.65 ± 0.15 mag 
color index. Broadly speaking, SED fitting is a trade-off between spectral type, luminosity 
and line-of-sight extinction. The assumption of identical extinction toward V 773 Tau C and 

V 773 Tau AB can be relaxed, in which case our estimated spectral type and luminosity are 
modified. With an additional lmag of extinction at V toward V 773 Tau C, we find that its 
SED is well fit by a K7 dwarf with a luminosity of 2.3 ± 0.3 L ; a significant excess at V 
remains present although the H and K band fluxes are then consistent with photospheric 
levels. Larger extinctions are not allowed, since it would imply that more than 0.8 L & is 
absorbed by dust and reemitted at longer wavelengths, while this is the total observed mid- 
to far-infrared luminosity in the system (see below). We adopt here the "equal- Ay" estimates 
for V 773 Tau C, i.e. a M0. 5 spectral type and a 1.5 L® luminosity. 

In contrast to these three components, which can be fitted with normal photospheric 
colors and a moderate NIR excess for V 773 Tau C, V 773 Tau D displays extremely red colors 
(H — K ~ 1.3 mag, K — L' ~ 1.7 mag). A blackbody fit to its dereddenned if, K, L' photo- 
metric data yields T NIR = 900 ± 200 K, which is much cooler than any stellar photosphere. 
Its 1.6-3.8 ^m SED is rising and this extends to longer wavelengths since V773 TauD is the 
brightest component of the system in our M' image, suggesting that it might the dominant 
source of the unresolved mid-infrared flux of the system. Assuming that the IRAS fluxes 
can be fully assigned to V 773 TauD, we conclude that the SED of this object presents a 
broad peak from ~ 5 fim to ~ 30 fim and that its bolometric luminosity is 0.8 ±0.1 L Q . Since 

V 773 Tau C may exhibit a significant excess beyond the M' band, either through reprocess- 
ing of absorbed visible starlight or through emission of its disk, this value is an upper limit to 
the luminosity of V 773 TauD. Nonetheless, we consider that the extremely red 1.6-4.7/zm 
colors of V 773 Tau D are substantial evidence that the mid-IR emission arises predominantly 
from this component. We also find that the spectrum of V 773 TauD is remarkable in that 
it shows no photospheric features, in spite of the high signal-to-noise ratio of our data (Fig- 
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ure3). The only significant feature in this spectrum is the Br7 line, which displays a weak 
emission with an equivalent width of 0.4 ± 0.1 A. Since all photospheres have identifiable 
features in the observed wavelength range, we conclude that the NIR light we are receiving 
from V773TauD has been thermally reprocessed by dust. 

Finally, we note that at least some of the components in this system are photometrically 
variable, potentially increasing the uncertainties on the luminosities derived above. For 
instance, the flux ratio between V773TauAB and V773TauC varies by 1-1.5 mag both 
in the visible and the NIR (see Table 1). In fact, the resolved photometry presented by 
Ghez et al. (1997) suggests that V773TauC is the variable component, as its K band 
magnitude varied by ~ 1 mag while V 773 Tau AB remained almost unchanged. Furthermore, 
V773 Tau AB, which dominates the flux of the system in the visible, has been shown to vary 
by only ~ 0.15 mag in the visible (Rydgren & Vrba 1983). Similarly, the V773TauAB to 

V 773 Tau D flux ratio at K band varied by ~ 1 mag between 1998 and 2002 implying that 

V 773 Tau D is also variable. In the absence of simultaneous visible and NIR photometry 
for all components of the system, it is not possible to decide whether the variability of 

V 773 Tau C and V 773 Tau D is due to varying extinction along the line of sight or is intrinsic 
to these objects. 



4. Discussion 
4.1. Nature of the components of the system 

4.1.1. V773TauAB and V773TauC 

V773TauAB has so far been studied spectroscopically at optical wavelengths, under 
the assumption that the contamination by V 773 Tau C is negligible. The weakness of the Ha 
emission line in its spectrum (2-4 A equivalent width; Herbig, Vrba & Rydgren 1986; Cabrit 
et al. 1990; Feigelson et al. 1994) as well as the chromospheric-like ultraviolet Mgll h and k 
emission lines (Feigelson et al. 1994), have led to a WTTS classification of both components 
of the SB. The observed Br7 line in absorption, the absence of NIR broadband excess in 
its SED up to L', and the very limited amount of visible and NIR variability (§3.2) are 
also consistent with this classification, and we therefore adopt it. Based on their estimated 
spectral types and luminosities, White & Ghez (2001) concluded that M A = 1.46 M Q . Our 
adopted luminosity is marginally larger, and we adopt M A ps 1.5 ± 0.1 M Q . Using the 
dynamical mass ratio estimated by Welty (1995) for the SB, this implies a mass of ps 1.1 M Q 
for V 773 Tau B. 
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The classification of V773TauC as a WTTS by White & Ghez (2001) was based on 
its photopheric-like SED. However, as pointed out in §3.2, it displays conclusive excesses at 

V and M' as well as possible excesses at H and K. Its NIR excess is typical of CTTS and 
unheard of in WTTS. Similarly, the observed variability of V773TauC in the visible and 
the NIR is typical of CTTS: variations of more than ~ 0.5 mag in the visible and/or the NIR 
are generally associated with actively accreting TTS (Herbst et al. 1994; Skrutskie et al. 
1996). It thus appears that V773TauC is a CTTS, and such a classification resolves some 
of the puzzling properties of this system. These include the broad Ha emission line found by 
Feigelson et al. (1994), the strong Br7 emission line 6 measured by Folha & Emerson (2001) 
and the detection by Carbit et al. (1990) of visible forbidden emission lines, which are all 
atypical for a WTTS and rather suggest CTTS-like ongoing accretion and mass loss. All of 
these features can be explained if V 773 Tau C is a normal CTTS whose spectrum has been 
diluted by the spectrum of its bright WTTS neighbour since it accounts for only 20-30 % 
of the flux of the unresolved system. We therefore conclude that V 773 Tau C is a CTTS, 
which should be confirmed in the future by obtaining separate spectra for V773TauAB and 

V 773 Tau C in the visible or the NIR. 

Estimating the mass of V 773 Tau C is a non-trivial task, as it lies above the 1 Myr 
isochrone of the evolutionary model of Baraffe et al. (1998), with either M0. 5 or K7 adopted 
as its spectral type, while V773TauA falls on the 5 Myr isochrone (White & Ghez 2001). 
This is most likely because the luminosity of V 773 Tau C is partially contaminated by 
accretion-induced emission and is not purely photospheric. White & Ghez estimated an 
M0 spectral type and M c = 0.7 M because the object lies in the HR diagram just above 
this almost vertical evolutionary track. Applying a similar reasoning, our K7-M0.5 esti- 
mated spectral type implies a stellar mass in the range 0.65-0.8 M & . We therefore adopt 
Mc = 0.7 ± 0.1 M Q . The total mass of the V773TauABC triple system is therefore pre- 
dicted to be 3.3 ± 0.2 M Q , which is in agreement with the dynamical mass we derived in 
§3.1. While a finer analysis of the stellar masses in the system will require radial velocity 
measurements, this indicates that the stellar masses we adopt here are coherent with our 
dynamical analysis. 



6 In a 1.3 /an spectrum obtained two nights earlier, they failed to detect the Pa/? line in emission, even 
though it is usually stronger than the Btj in CTTS (Muzerolle et al. 1998). This probably indicates that 
the emission lines are strongly variable. 
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4-1-2. V773TauDisan "infrared companion" 

As noted in §3.2, the SED of V773TauD cannot be fit with any stellar photosphere 
and its spectrum suggests that this object is seen through optically thick dust. Therefore, it 
cannot be classified as a usual TTS. In fact, the shape of its SED matches that of a handful 
of objects first identified as IRCs to TTS by Chelli et al. (1988) and extensively studied by 
Koresko, Herbst & Leinert (1997). These companions were initially defined as extremely red, 
cool, objects located within a few hundred AU of a TTS. In agreement with our observations 
of V773TauD, these IRCs have SEDs that are too broad to be fitted by single blackbody 
curves and that dominate the flux from the systems at wavelengths longer than 2-5 /im. 
Furthermore, K band spectroscopy of several other IRCs have revealed featureless spectra, 
with the exception of Br7 in emission (Beck, Prato & Simon 2001; Duchene et al. 2002; 
Herbst, Korseko & Leinert 1995). Overall, we consider that the available dataset regarding 

V 773 Tau D clearly supports its classification as a new member of the small category of IRCs. 

Although V 773 Tau D shares enough properties with other IRCs to be classified as 
such, it stands out from the other members of this class in at least two respects. The first 
remarkable property of V 773 Tau D is that its bolometric luminosity is not larger than that 
of the system's primary. In fact, all other components in the system are intrinsically brighter 
than the IRC. The upper limit on its bolometric luminosity places a strict upper limit on its 
mass of about 1.2 M (Baraffe et al. 1998, Palla & Stahler 1999, Siess et al. 2000), possibly 
significantly less if its luminosity is partially driven by accretion. Following the analysis 
of Do-Ar24E by Koresko et al. (1997), we believe that the low luminosity of V 773 Tau D 
indicates that its mass is smaller than that of the other components in the system. We 
therefore conclude that M D < 0.7 M Q . 

The other noticeable property of V 773 Tau D can be found in the 8-13 /iin spectrum of 
the unresolved system presented by Hanner, Brooke & Tokunaga (1998), which displays no 
silicate feature, either in absorption or in emission. If the mid-infrared flux of V 773 Tau is 
dominated by the IRC, this contrasts with all other IRCs for which a spectrum has been 
obtained (Hanner et al. 1998; Natta, Meyer & Beckwith 2000) since they all present a strong, 
broad, spectral feature around 10 /im. The absence of a silicate feature in the spectrum of 

V 773 Tau D is unlikely to be the result of a strong silicate depletion in its envelope, since this 
has not been observed among TTS. Possible explanations for this phenomenon include time 
variability of the feature, such as observed in the case of another IRC (Glass I, Giirtler et al. 
1999), an accidental though almost exact cancellation of emission and absorption features in 
the envelope, the growth of most silicate grains beyond a few microns in size, or an unusually 
low temperature in the outer parts of the envelope. In cold environments (T < 100 K), water 
ice can condense on and shield small silicates grains, in which case the grain opacity in the 



-14- 



mid-infrared does not show the usual marked features at 10 /mi and 20 /tin (e.g., Pollack et 
al. 1994). It is reasonable to expect that the outer parts of the dusty envelope surrounding 
V 773 Tau D are colder than that of other IRCs as it is a less luminous source, though the 
actual cause of the atypical mid-infrared spectrum of this source remains an open issue. 

4.1.3. Nature of the IRCs 

The exact nature of IRCs is still debated. They may be normal TTS undergoing episodes 
of intense accretion and which are embedded in optically thick dusty envelopes with radii as 
small as a few AU (Ghez et al. 1991; Koresko et al. 1997; White & Ghez 2001; Duchene et al. 
2002). They also may be embedded protostars, i.e. objects in an earlier evolutionary phase 
(e.g., Ressler & Barsony 2001). This latter interpretation would imply that IRC binary sys- 
tems are not coeval, in contradiction with results on large samples of TTS binaries (Hartigan, 
Strom & Strom 1994; White & Ghez 2001). We emphasize that the observational distinction 
between these two types of objects is quite limited; both are highly embedded sources whose 
flux is dominated by accretion onto and/or contraction of the central object (White & Ghez 
2001), from which significant photometric variability can be expected. Photometric variabil- 
ity has indeed been documented for several IRCs (Ghez et al. 1991, Koresko et al. 1997), 
including V 773 Tau D (§3.2). Because IRCs are, by definition, encountered near few Myr- 
old TTauri stars, it is difficult to imagine that they recently formed from the contraction 
of a vast (several thousand AU) envelope as normal TTauri stars do, casting some doubt 
upon their protostellar nature. However, an isolated object with the same photometric and 
spectroscopic properties as V 773 Tau D would unambiguously be classified as an embedded 
protostar. The future detection of photospheric features in the spectrum of V 773 Tau D 
might clarify its properties, but for now we cannot conclusively determine the nature of this 
object. 

Among the seven IRCs studied by Koresko et al. (1997), six were found to be more mas- 
sive than their optically bright "primaries", which suggested that their peculiarly dusty envi- 
ronment were a consequence of being the most massive component in the systems. However, 
the new IRC discovered here within V 773 Tau is found to be the lowest luminosity /mass 
object of the quadruple system. We therefore propose that there is a continuum of IRC 
masses and luminosities. In other words, there may be many yet undetected IRCs, as obser- 
vational studies of such objects have likely been biased towards the brightest and therefore 
most massive ones. If so, we may expect that further deep, high-angular resolution 10-20 /im 
imaging around known TTS will reveal a significant population of much fainter objects with 
SEDs similar to that of V 773 Tau D. 
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4.2. A "fully mixed" system 

Most IRCs studied so far are associated with CTTS (Koresko et al. 1997). However, 
the quadruple system V773Tau consists of two 1-1.5 M & WTTS in a SB, a 0.7 M CTTS 
orbiting around it at ~ 15-25 AU, and a low-mass IRC located further away. This system can 
therefore be considered as a "fully mixed" system, containing all three types of TTS known 
so far: CTTS, WTTS and IRC. How it has evolved in the last Myr or so remains mysterious 
and it may even seem puzzling that a multiple system with two WTTS can also contain 
more active objects such as a CTTS and an IRC. However, one can speculate that strong 
interactions between the various components and their circumstellar environment have been 
major factors in its evolution. 

The presence of two close companions located ~0.3 AU and ~15 AU away from V 773 Tau A 
can naturally explain why the two components of the SB have rapidly evolved toward WTTS. 
The various orbital motions would indeed clear out any circumstellar disk except in a nar- 
row annulus between ~ 1-5 AU (e.g., Artymowicz & Lubow 1994) within a few thousand 
years. Similarly, the circumstellar disk surrounding V 773 Tau C would be rapidly truncated 
to an outer radius of at most 3-4 AU. While the material surrounding V773TauAB appears 
to have fully dissipated, V 773 Tau C still possesses an active accretion disk. However, the 
lifetime of such a small disk is relatively short unless it supports only a very small accretion 
rate. The upper limit on the millimeter thermal emission from the system corresponds to a 
maximum possible disk mass of only ~ 10~ 2 M (Beckwith et al. 1990; Dutrey et al. 1996). 
Given the age of the system and assuming that all this mass is located around V 773 Tau C, 
only a very low time-averaged accretion rate (< M > < 2 x 10~ 9 M /yr) could have allowed 
the disk to survive for so long. Most likely, this circumstellar disk is being replenished from 
an outer reservoir of material, as has already been suggested for tight pre-main sequence 
binaries (e.g., Prato & Simon 1997). The fact that V 773 Tau C seems to be the only com- 
ponent of the central triple system to be replenished may suggest that the orbital motion 
of the SB prevents material from falling deep into the potential well of V 773 Tau AB due to 
resonance interactions with the two stars. 

Since it is on a much wider orbit, V 773 Tau D can retain a much larger envelope, several 
tens of AU in diameter, and its dissipation timescale may not be shorter than the age of the 
system. Although replenishment is not required to account for the presence of circumstellar 
material around this source, it is nonetheless surrounded by an unusually optically thick 
envelope. This may be a transient state generated through interactions with the other 
components of the system. In any case, V 773 Tau is the first system among TTauri stars 
to display such a range of evolutionary states and it is likely that the compactness of the 
system, through strong star-star interactions, has an important role to play in this variety. 
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4.3. Do passive disks around TTS really exist? 

The 5-60 /im flux from V 773 Tau, that was previously assigned to a truncated, fossil, 
circumstellar/circumbinary disk around the WTTS system V773TauAB, is now believed to 
be associated with the newly identified IRC, V 773 Tau D, with a possible contribution of the 
disk surrounding the CTTS V 773 Tau C. We conclude that there is no evidence for a fossil 
disk in this system. 

Candidate fossil disks among pre-main sequence stars are quite rare. Skrutskie et al. 
(1990) and more recently Stassun et al. (2001) obtained 10 fxm photometry for relatively large 
samples (83 TTS and 32 WTTS respectively) and only identified a handful of such cases 
among pre-main sequence objects in Taurus and Orion. The most suggestive candidate is 
GMAur which is known to be surrounded by a dusty disk in Keplerian rotation (Koerner, 
Sargent & Beckwith 1993) but shows at best a marginal photometric excess shortwards of 
10 /im (Stassun et al. 2001, White & Ghez 2001). However, this object presents a strong 
accretion-induced ultraviolet /blue excess and it is unclear whether it should be classified as 
a passive disk. About ten other objects, or roughly 10 % of the combined sample, might fall 
in this category, suggesting a lifetime for this phenomenon on order of a few 10 5 yrs. 

However, our analysis of V 773 Tau has led us to remove this system from the list of 
candidate fossil disks and one may wonder whether other systems will similarly be reclassified 
out of this category by further analysis. If so, fossil disks would be extremely short-lived 
structures. As we have shown here, the presence of a close, unresolved IRC in the vicinity 
of a WTTS can be misinterpreted as a passive disk on the basis of its SED. Similarly, a 
tight SB would rapidly clear a central gap in its surrounding disk and the absence of the 
hottest parts of the disk would result in a dip in the NIR part of the SED of the system 
(e.g., Mathieu et al. 1991). A disk inner radius of about 0.4 AU results in no excess above 
photospheric levels up to 5 fim or so, similar to the SED of a passive disk. Such an inner 
radius can be created by a binary with a semi-major axis of about 0.15 AU (Artymowicz 
& Lubow 1994), corresponding to an orbital period of about 20c? for a 1 M & total system 
mass. Unresolved IRCs and short period SBs can therefore mimic passive disks, and it is 
therefore important to compare the frequencies of these various systems to decide whether 
a significant proportion of passive disks really exists among pre-main sequence stars. 

So far, at least five IRCs out of ~ 100 pre-main sequence objects are known in Taurus 
(TTauSa, Haro6-10N, XZTauB, UYAurB and now V 773 Tau D) and more may be dis- 
covered when 5-10 /im high-angular images of large samples of TTS are performed. This 
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observed proportion of ~ 5 % of IRCs is therefore only a lower limit 7 . The proportion of 
SBs with P orb < 20 d for solar-type main sequence objects is estimated to be about 4-5 % 
(Duquennoy & Mayor 1991), and it might be even higher among young stars in Taurus as is 
the case for visual binaries (e.g., Ghez et al. 1993). To the best of our knowledge, only 4 such 
short-period SBs are known among the Taurus pre-main sequence population, out of which 
only 2 (V826Tau, Mundt et al. 1983, and LkCa3, Mathieu 1994) are classified as WTTS 
based on their visible/NIR photometric and spectral properties. Overall, the proportion of 
passive disks (< 10 %) qualitatively matches those of IRCs and short period SBs among 
pre-main sequence stars (both > 5%). It is very possible that most (if not all) fossil disks 
candidates are in fact associated with currently unknown IRCs or SBs. More generally, this 
emphasizes how dangerous SED analyses are if all stellar components have not been resolved. 
This issue, which was suggested by Chelli et al. (1988) and Ghez (1996), resurfaces as a 
concern with the advent of IR missions with only moderate spatial resolution (> 1"), such 
as SIRTF and SOFIA. 



5. Conclusion 

Using speckle interferometry techniques at the 10m Keck I and 200-inch Hale telescopes 
and the AO system on Keck II, we have obtained multi-epoch high-angular resolution images 
of the multiple system V 773 Tau in the visible and the NIR throughout the last decade. It 
had been considered one of the rare examples of the so-called "passive" or "fossil" disks 
which are thought of as disks in which planetesimals have formed in the inner parts, thus 
halting the accretion process on the stars without expelling the outer disk material. 

In addition to the unresolved SB V 773 Tau AB and the previously known speckle com- 
panion V 773 Tau C, we identify a K ~ 9, K — L' ~ 1.7 mag, fourth component in our 
AO images and in some speckle datasets. V 773 Tau D, which is located only 0'.'2 away 
from V773TauAB, is the brightest source in the system at 4.8 fxm. We have also obtained 
medium-resolution 2 fxm spectra of V 773 Tau AB and V 773 Tau D. Given its NIR brightness, 
its location, and its apparent motion, we conclude that the newly identified component is 
physically associated with the system. We fit the motion of V 773 Tau C around V 773 Tau AB 
with a closed Keplerian orbit, and we conclude that the quadruple system can be dynami- 
cally stable provided V773 TauD is on a highly inclined orbit, as are the other orbits within 



7 IRCs constitute a significant fraction of the "high accretion rate" population identified by White & Ghez 
(2001) in Taurus on the basis of the extremely red colors of these objects (K — L > 1.4, as is the case for 
V 773 TauD). Roughly 10% of all TTS belong to this other category, in qualitative agreement with our 
estimated proportion of IRCs. 
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the system. As opposed to TTau (Loinard et al. 2003), this system does not support the 
proposed explanation to the observed overabundance of young high-order multiple systems 
that many of them are in fact unstable systems that will rapidly decay into an ejected single 
star and a lower order multiple system (Ghez et al. 1993; Reipurth 2000). 

Analysis of the SED of the system and of our spectra supports the following picture: 
both components of V 773 Tau AB are 1-1.5 M WTTS, V 773 Tau C is a lower-mass variable 
CTTS and V 773 Tau D is a new member of the IRC category, i.e. a deeply embedded 
companion to a normal TTauri star. The properties of this new object include a broadly 
peaked SED from ~ 5-20 fim, significant variability at 2 /im, and a featureless 2 fim spectrum. 
Its low bolometric luminosity (Lp < 0.8 L©) shows that it is the least massive component 
of the system, in contrast to most other IRCs. We suggest that many low-mass IRCs might 
remain to be discovered. The physical nature of IRCs remains unconstrained since both 
embedded TTauri stars or less evolved protostars are almost identical in terms of their 
observable properties. 

The immediate implication of the presence of V 773 Tau D is the rejection of the hy- 
pothesis of a passive disk in the system since the long known IR excess is now attributed 
to the IRC. Given the estimated proportions of IRCs and short period SBs, we argue that 
the frequency of passive disks among pre-main sequence stars may have been strongly over- 
estimated in the past and that the lifetime of this transition phase may be much shorter 
than previsouly thought. More generally, our analysis of V 773 Tau shows how cautious one 
must be when using NIR or mid-IR photometry of TTS obtained with low spatial resolution 
devices. 
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Table 1. Summary of all spatially resolved measurements of the V773Tau system. 



Date 


Obs. 


Filter 




V773TauAB- 


-C 




V 773 Tau AB- 


-D 






Tech. 




p (mas) 


9 (°) 


Am (mag) 


p (mas) 


6 (°) 


Am (mag) 


Ref. 


1990 Oct 03 


Spk 


K 


112±2 


295±4 


0.80+0.03 


T—s 

"' 




T-^: 


2 


1991 Sep 21 


lD a 


K 


170±10 


295±3 


2.22+0.33 




(") 


( ) 


3 


1991 Sep 25-29 


Spk 


-^620 


106±5 


298±2 


2.94+0.10 


(") 
(...) 


(") 
(...) 


( ") 
(...) 


1 


55 




^656 


55 


55 


2.19+0.10 


(...) 


(...) 


(...) 


1 


55 




^800 


55 


55 


2.05+0.10 


(...) 


(...) 


(...) 


1 


55 




^850 


55 


55 


2 19+0 10 


(...) 


(...) 


(...) 


1 


1992 Oct 10 


Spk 


K 


120±20 


307±5 


1.71+0.01 


(...) 


(...) 


(...) 


4,5 


1993 Oct 05 


Spk b 


J 


111+4 


308 9+1 4 


2 37+0 05 


(...) 


(...) 


(...) 


6 


1993 Nov 25 


Spk 


K 


92±8 


304±8 


1.75+0.22 








4,5 


1994 Oct 19 


Spk 


K 


65.6±1.4 


318±2 


0.46+0.04 


(") 


(") 


,"\ 


4 5 

^5 u 


1994 Oct 29 


HST 


jji 


62.8±2.4 


321 9+2 7 


1 S9+0 78 


(") 
(...) 


(") 
(...) 


('") 
(...) 


5 


55 


55 


B' 


55 


57 


2.78+0.92 


(...) 


(...) 


(...) 


5 


55 


55 


V 


55 


11 


i VQ+n 44 








5 


55 


55 


R! 


55 


11 


1.36+0.37 








5 






r 




11 


1.19+0.54 








5 


1996 Oct 24 


Spk c 


K 


(...) 


(...) 


(...) 


162+29 


176+5 


2.17+0.27 


1 


1997 Nov 19 


Spk c 


K 


(...) 


(...) 


(...) 


187.7+3.9 


173.7+1.5 


2.12+0.10 


1 


1998 Nov 04 


Spk 


K 


57.4±1.2 


67.1+1.4 


0.76+0.03 


191.3+3.6 


164.7+1.3 


2.85+0.06 


1 


1999 Nov 20 


Spk 


K 


74.4±1.4 


81.9+1.1 


0.83+0.06 


202.4+3.8 


161.3+2.0 


2.79+0.05 


1 


2000 Nov 20 


AO 


H 


85±3 


88.4+1.5 


0.44+0.10 


212+4 


155.7+0.7 


2.86+0.15 


1 






K 






0.46+0.05 






1.92+0.10 


1 


2000 Dec 08 


Spk 


K 


91.1±3.9 


88.5+1.2 


1.01+0.07 


217.4+4.9 


159.3+4.9 


2.43+0.07 


1 


2001 Feb 08 


Spk d 


K 


95±5 


93.5+0.6 


(...) 


(...) 


(...) 


(...) 


7 


2001 Nov 02 


Spk d 


K 


99±7 


94.6+1.0 


(...) 


(...) 


(...) 


(...) 


7 


2001 Dec 08 


Spk 


K 


100.0+1.0 


92.3+1.0 


0.33+0.07 


223.0+3.0 


154.1+1.0 


1.84+0.07 


1 


2002 Dec 13 


AO 


L> 


116.5±2.0 


100.6+0.7 


0.06+0.05 


237+3 


153.7+0.6 


0.36+0.05 


1 


55 


55 


M' 


55 


11 


-0.21+0.05 


55 


11 


-0.30+0.05 


1 



Note. — The data presented in this table have been obtained through two-dimensional ("Spk") and one-dimensional ("ID") 
speckle intcrferometry, AO or i?5T/WFPC2 imaging, p represents the binary separation and its position angle as measured 
East from North on the sky. In all observations between 1990 and 1997, only one companion was unambiguously detected. 

a This measurement was obtained by analyzing two orthogonal one-dimension scans and retrieving the actual separation and 
position angle by deprojection from these axes. The uncertainties associated with this method are significantly larger from 
those associated with two-dimension speckle measurements, and this data point differs by more than 5<r from the visible 2D 
speckle result we obtained just a week later, suggesting that the uncertainties of the ID measurement are underestimated. It 
has not been included in our final analysis. 

b Instead of being a Nyquist-samplcd measurement, this dataset was obtained with a pixel scale that is almost equal to the 
diffraction limit of the 3.5 m Calar Alto telescope in the J band (0"07/pixcl and 0'.'074 respectively). This measurement is thus 
significantly aliased and therefore has not been used in our analysis. 

c In these measurements, the tight pair V 773 Tau AB— C was not fully resolved. All photometric and astrometric information 
extracted from this dataset therefore pertains to the V 773 Tau ABC-D system. 

d These measurements were obtained with a different instrumental and calibration set-up than most of our data. To minimize 
the effect of systematic differences between datasets, we do not include these datascts in our analysis since our own data already 
provide an adequate time sampling throughout the 1998-2002 period. 



References. — 1) This work; 2) Ghez et al. (1993); 3) Lcincrt ct al. (1993); 4) Ghcz ct al. (1995); 5) Ghcz ct al. (1997); 6) 
Woitas ct al. (2001); 7) Tamazian et al. (2002). 



Table 2. Parameters for the orbit of V 773 Tau C around V 773 Tau AB. 



M ABC (M Q ) P orfe (yrs) a (AU) e T u (°) 

3.7±0.7 a 46 ±6 20 ±1 0.3 ±0.1 66 ± 3 b 1996.5 ± 0.8 81 ± 10 b 288 ± l b 



a The uncertainty on the total system mass only includes the fitting uncertainty. An 
additional 10 % random uncertainty arises from the error on the distance to the system since 
orbital analysis relying on astrometric data only constrains M syst /D 3 . 

b In the absence of radial velocity measurements, an ambiguity in the inclination angle re- 
mains unsolvable. We have assumed that % < 90°, which means that V 773 Tau C is currently 
in front of the other component. However, it is equally possible that i = 114°, in which case 
iv = 279° and Q = 108°. The uncertainties on these parameters, as well as the location of 
periastron, are however unaffected by this ambiguity. 
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1 arc se c 

Fig. 1. — Examples of speckle interferometry power spectra of the V 773 Tau system obtained 
in 1997 (a), 1999 (b) and 2001 (c). In the first epoch, only the V 773 Tau AB-D pair, oriented 
more or less North-South is clearly detected. In the later two epochs, a set of much wider 
fringes appears, corresponding to a tight pair oriented roughly East-West, V773TauAB- 
C. In the 2001 dataset, it is particularly clear that two independent sets of fringes are 
detected or, equivalently, that the system is resolved into three different components. Note 
the different spatial scales for the Palomar (a and b) and Keck (c) datasets, which explains 
why the fringes associated to the V773TauAB-C pair in the latter show several cycles in 
the Fourier space instead of merely one. 
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V 773 Tau C 



Dec. 2002 



Fig. 2. — Adaptive optics images of the V 773 Tau multiple system obtained in November 
2000 (top row) and December 2002 (bottom row). A square root stretch has been used in 
all images. North is up and East to the left and each image is 1" across. The lowest contour 
in the H and K images is about 120 times and 150 times higher than the background noise. 
The contours represent 2 %, 3 %, 5 %, 15 %, 33 %, 60 % and 90 % of the peak value in the H 
band image and 3%, 7%, 15 %, 45 % and 90 % of the peak at K. The motion of V 773 Tau C 
with respect to V 773 Tau AB between the two epochs is particularly obvious in these images, 
as well as the relative brightening of V 773 Tau D towards longer wavelengths. 
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Fig. 3. — Continuum-normalized spectra of V773TauAB and V773TauD; the latter has 
been shifted upwards by 0.05 for clarity. The arrows indicate all known spectral features in 
this wavelength range (Kleinmann & Hall 1986) and the Br7 line is indicated by the longer 
arrow. Both spectra have been smoothed down to a resolution of about R ~ 1900. A few 
spikes, both in emission and in absorption, appear in the spectrum of V 773 Tau D but they 
are too narrow (typically 1 pixel-wide) to be real features. Furthermore, similar "features" 
can be seen in the spectrum of V773TauAB (e.g., at 2.217/im or 2.267 /xm), which suggests 
that they correspond to poorly corrected telluric absorption lines. 
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Fig. 4. — Location of the companions that have been spatially resolved from V 773 Tau AB 
since 1990. The SB (large cross) is fixed at location (0,0); the amplitude of the internal 
orbital motion of V 773 Tau AB is negligible at the scale of this figure. The length of the 
tick marks indicate 1 a error bars and the year of observation is associated to each point. 
Only the datapoints used in the dynamical analysis are shown here. Our best fit orbit for 
V 773 Tau C is represented by the solid curve, together with its periastron (open diamond) 
and line of nodes (dashed line). 
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Fig. 5. — Dereddened SED of the components of the system. Diamonds: unresolved system; 
triangles: V773TauAB; squares: V773TauC; asterisks: V773TauD. The dashed curve is 
the sum of the K2 and K7 photospheres fitted to the components of V 773 Tau AB while the 
dotted curve is the MO. 5 photosphere fit to V 773 Tau C. The dot-dashed curves joins the ob- 
served HKL' photometric points for V 773 Tau D with the unresolved IRAS measurements. 



